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ABSTRACT 



The GOV dwarf TrES-2A, which is transited by a hot Jupiter, is one of the main short-cadence targets of the Kepler telescope and, 
therefore, among the photometrically best-studied planetary systems known today. Given the near-grazing geometry of the planetary 
orbit, TrES-2 offers an outstanding opportunity to search for changes in its orbital geometry. Our study focuses on the secular change 
in orbital inclination reported in previous studies. We present a joint analysis of the first four quarters of Kepler photometry together 
with the publicly available ground-based data obtained since the discovery of TrES-2b in 2006. We use a common approach based on 
the latest information regarding the visual companion of TrES-2A and stellar limb darkening to further refine the orbital parameters. 
We find that the Kepler observations rule out a secular inclination change of previously claimed order as well as variations of the 
transit timing, however, they also show slight indication for further variability in the inclination which remains marginally significant. 

Key words, stars: individual: TrES-2 - stars: planetary systems 



1. Introduction 

The transiting Jovian planet TrES-2b was originally discovered 
by O' Donovan et al. (2006), and later extensively observed by 
Holman et al. (2007) and Sozzetti et al. (2007), who determined 
more accurate properties for the TrES-2 system. TrES-2A is an 
old 5 Gyr) GOV star with about solar radius and mass, or- 
bited by a close-in (0.035 AU) planet with a mass of 1.2 My and 
a period of 2.47 d. The planet has a radius of 1 .24 R j and a large 
inclination angle of about 84°, implying a nearly grazing tran- 
sit. The spectroscopic parameters of the systen determined by 
Sozzetti et al. (2007) were confirmed and refined by Ammler- 
von Eiffetal. (2009). 

As already pointed out by Holman et al. (2007), TrES-2 of- 
fers an outstanding opportunity to search for changes in orbital 
parameters for several reasons: Its large inclination and close-in 
orbit yield a high impact parameter, making the shape and tim- 
ing of the transit sensitive to variations of the orbital parameters. 
Its large radius and short orbital period allow the study of large 
numbers of deep transits over hundreds of epochs. Additionally, 
TrES-2 is in the field of view of the Kepler telescope, which 
has already gathered more than 220 days of publicly available 
photometric data, probably making TrES-2 the photometrically 
best-studied exoplanetary system to date. 

Since its discovery there have been numerous attempts to de- 
tect secular changes in the orbital parameters of TrES-2b due to 
an external perturber. Mislis & Schmitt (2009) combined their 
own transit observations with the data from Holman et al. (2007) 
and found these data to be consistent with a decrease in transit 
duration, suggesting a decrease of about 0.075°/yr in the orbital 
inclination of the planet, an interpretation corroborated by an- 
other series of observations (Mislis et al. 2010). Reverting to sec- 
ular perturbation theory, Mislis et al. (2010) suggest that a Jovian 
mass planet with a period between 50 to 100 days can explain the 
observed inclination changes. On the other hand, Scuderi et al. 
(2010), in their analysis from data taken in June 2009, find sys- 



tematically higher values for the inclination of TrES-2b, ruling 
out any systematic trend. Their analysis also includes 33 days of 
Kepler data comprising multiple transits originally published by 
Gilliland et al. (2010), which they find to be in agreement with 
their ground-based observations. 

Additional photometric data of TrES-2 are presented by 
Rabus et al. (2009), who analyze five transits sytematically 
searching for signatures of transit timing variations. However, 
they do not find a clear signal and are only able to provide upper 
constraints on the perturber's mass. Christiansen et al. (2011) 
obtained nine additional transits from NASA's EPOXI Mission 
of Opportunity (Deep impact). Comparing their results with ear- 
lier measurements, they find a formally significant decrease in 
the transit duration with time, which however disappears if they 
exclude the duration measurement by Holman et al. (2007). 

A thorough analysis of the first two quarters of Kepler data 
(Q0 and Ql) is presented by Kipping & Bakos (201 lb), who do 
not find any inclination change within the Kepler data, obtain- 
ing an inclination change of (+0.0019 + 0.0020)°/cycle. They 
also reanalyze the data by Holman et al. (2007) fitting the limb- 
darkening coefficients and compare the measured average transit 
durations in both data sets. Their results neither support nor ex- 
clude the inclination change, largely due to the still very small 
temporal baseline of 18 cycles. Raetz et al. (2009, 201 1) analyze 
22 additional transits observed between March 2007 and August 
2010, but find no indications of variations in the transit dura- 
tion. However, their data quality does not suffice to exclude the 
inclination change proposed by Mislis et al. (2010). 

The whole discussion is further complicated by the results 
of Daemgen et al. (2009), who detected a visual companion to 
TrES-2A at a separation of 1.089 + 0.008" that contaminates 
the photometric data, and leads to filter-dependent deviations in 
the orbital parameters if not being accounted for. While Scuderi 
et al. (2010) already mention the contribution of a third light 
and estimate its influence on the system parameters, Kipping & 
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Bakos (201 lb) were actually the first to include it in their light- 
curve modeling of the Kepler data. 

The purpose of this work is to present a joint analysis of all 
publicly available Kepler data of TrES-2 focussing on the ap- 
parent change in orbital inclination as proposed by Mislis et al. 
(2010). We exploit four quarters of Kepler photometry and pub- 
licly available ground-based transit observations in a joint ap- 
proach including the latest observational results regarding the 
TrES-2-system, aiming at a consistent picture of the observa- 
tional status quo regarding changes in the orbital parameters, es- 
pecially the inclination of the planet. 

2. Light-curve modeling 

We fit the transit data with the transit model developed by 
Mandel & Agol (2002) making use of their occultquad and 
occultnl FORTRAN routines 1 . From the transit light curve, we 
can directly infer the following parameters: the orbital period, P, 
the time of transit minimum, t, the radius ratio, p = R p /R s , 
the semi-major axis in stellar radii, a/R s , the orbital inclina- 
tion, i, and the limb-darkening law. For our fits we either as- 
sume a quadratic limb darkening with coefficients u\ and u 2 or a 
four-parametric non-linear prescription as proposed by Claret & 
Bloemen (201 1, their Eq. 5), which provides a significantly more 
accurate representation of model atmophere intensities (Howarth 
2011a). Furthermore, our model includes the contribution of a 
third light, L3, given as a fraction of the host-star's luminosity. 

We assume a circular orbit fixing e — 0, as there are no indi- 
cations for any non-zero eccentricity (Kipping & Bakos 201 lb). 

2.1. On the role of limb darkening 

There is an ongoing debate on the role of limb darkening in the 
analysis of planetary transit light curves regarding the question 
whether the limb-darkening coefficients (LDCs) should be fit- 
ted together with the remaining parameters or left fixed during 
the fit. On the one hand, LDCs are the only quantities in the pa- 
rameter set describing the transit that can be predicted theoreti- 
cally, given some basic spectroscopic information of the object. 
As such they are inherently dependent on the quality of available 
spectroscopic data and the subtleties of underlying assumptions 
entering the model calculations itself. Fitting the LDCs therefore 
assures that the resulting transit parameters are independent of 
the spectroscopic background model and that the uncertainties of 
the spectroscopic parameters are readily included into the transit 
model. Furthermore, the inclusion of LDCs as free parameters 
in the fitting process leads to considerably larger uncertainties 
in the remaining model parameters, reflecting the fact, that the 
LDCs are highly correlated to most parameters. 

On the other hand, the data quality of most ground-based 
measurements is not sufficient to substantially constrain the 
LDCs. Several fit parameters, especially the LDCs, are known 
to be mutually correlated, as demonstrated analytically by Pal 
(2008) and numerically by Southworth (2008) in the specific 
case of TrES-2. The mutual correlations among the parameters 
are particularly strong in the case of TrES-2, where the transit is 
nearly grazing, because TrES-2b only covers the outermost parts 
of the stellar disk that are dominated by strong limb-darkening 
effects. The remaining orbital parameters thus depend crucially 
on the details of the adopted limb darkening resulting in strong 
correlations with the LDCs. The amount of information on the 
limb darkening, which can be extracted even from the Kepler 

1 http : //www . astro . Washington . edu/user s/agol 



Table 1. Spectroscopic parameters of TrES-2 





Sozzetti et al. (2007) 


Ammler-von EifFet al. (2009) 


Teff (K) 


5850 + 50 


5795 ± 73 


logg (cgs) 


4.4 ±0.1 


4.30 ±0.13 


St (km/s) 


1.00 ±0.05 


0.79 ±0.12 


[M/H] 


-0.15 ±0.10 


0.06 ± 0.08 



data, is limited by the fact, that the transit of TrES-2b is nearly 
grazing, i.e. /i « 1. Assuming a quadratic limb-darkening law 
this approximation yields 

= 1 - («i + u 2 ) + p(m + 2m 2 ) - 0(p 2 ). (1) 

Therefore, to zeroth order we are only sensitive to the sum of the 
LDCs, i.e. the amount of deducible knowledge on the individ- 
ual parameter is inherently limited. This introduces a degener- 
acy between both LDCs due to their strong correlation. While, in 
principle, it would always be preferable to fit for LDCs, thereby 
cross-checking spectroscopic results and theoretical modeling, 
analyses of ground-based observations of TrES-2 cannot afford 
to refrain from including as much independently obtained infor- 
mation as possible. 

Therefore, in a first attempt, we will fit the Kepler data 
without assuming any a priori information on the LDCs (see 
Sect. 4.1) and compare our results to theoretical predictions. 
However, comparing the results from the Kepler data to several 
ground-based transit observations requires a common treatment 
of the data. We therefore eventually use fixed LDCs for all data 
sets assuming the theoretically predicted values (see Sect 4.2). 

The spectroscopic studies by Sozzetti et al. (2007) and 
Ammler-von Eiff et al. (2009) concordantly suggest that 
TrES-2A has an effective temperature of 5800 K, a surface grav- 
ity of 4.4, mictroturbulence velocity of 1 km/s and solar metallic- 
ity (see Table 1). Claret & Bloemen (201 l)provide LDCs in sev- 
eral filter bands calculated for a grid of PHOENIX and ATLAS9 
model atmospheres, employing two different calculation meth- 
ods. We linearly interpolated on the data grid to obtain a set 
of LDCs pertaining to a PHOENIX model closely matching the 
inferred spectroscopic properties of TrES-2A. Our interpolated 
LDCs for the different filter bands are summarized in Table 2. 

2.2. On the contribution of a third light 

Daemgen et al. (2009) report the detection of a faint visual com- 
panion of TrES-2A in their ground-based high-resolution im- 
ages obtained with the AstraLux lucky imaging camera at Calar 
Alto. The angular distance from TrES-2A to this object is only 
1.089 ± 0.009". It contaminates all previous published observa- 
tions and introduces a systematic error in the transit parameters, 
especially the transit depth. According to Southworth (2010), 
roughly 5 % of third light can be compensated for by decreas- 
ing a/R s by 1 %, p by 3 % and ; by 0.1 degrees. 

Daemgen et al. (2009) measured magnitude differences in 
the Sloan Digital Sky Survey (SDSS) z' and i' filters, that can be 
used to obtain the flux contribution of the third light, L3. They 
find Az' = 3.429 + 0.010 and Ai' = 3.661 + 0.016 resulting in 
flux contributions of 0.0425 + 0.0004 and 0.0343 + 0.0005 in the 
z' and i' bands, respectively. We assume that the nearby com- 
panion has an effective temperature of 4400 K, surface gravity 
of 4.5, solar metallicity and microturbulent velocity of 2 km/s 
compatible with its spectral type being K4.5 to K6 as derived by 
Daemgen et al. (2009). We used Gray's SPECTRUM together 
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Table 2. Non-linear limb-darkening coefficients obtained by lin- 
ear interpolation on the tables provided by Claret & Bloemen 
(201 1) based on PHOENIX atmosphere models. The last column 
gives the fraction of flux due to the visual companion detected 
by Daemgen et al. (2009). 



All MCMC calculations make extensive use of routines of 
PyAstronomy 4 , a collection of Python routines providing a con- 
venient interface to fitting and sampling algorithms implemented 
in the PyMC (Patil et al. 2010) and SciPy (Jones et al. 2001) 
packages. 



Filter 


ay 


a 2 


a 3 


U4 


L 3 


Kepler 


0.473° 


0.192° 






0.0276 


Kepler 


0.801 


-0.676 


1.053 


-0.388 


0.0276 


Stromgren u 


0.910 


-1.006 


1.388 


-0.355 


0.0045 


Stromgren v 


0.863 


-0.767 


1.872 


-0.617 


0.0102 


Stromgren b 


0.698 


-1.206 


1.572 


-0.599 


0.0180 


Stromgren y 


0.758 


-0.663 


1.203 


-0.454 


0.0227 


Johnson U 


1.085 


-1.412 


1.721 


-0.466 


0.0053 


Johnson B 


0.754 


-0.937 


1.734 


-0.641 


0.0137 


Johnson V 


0.742 


-0.602 


1.134 


-0.434 


0.0216 


Cousins R 


0.796 


-0.632 


0.993 


-0.378 


0.0294 


Cousins I 


0.853 


-0.746 


0.963 


-0.363 


0.0360 


SDSS i' 


0.829 


-0.688 


0.935 


-0.354 


0.0343* 


SDSS z' 


0.861 


-0.807 


0.967 


-0.360 


0.0425 fc 



Notes. ( °' Interpolated coefficients for a quadratic limb-darkening law. 
(i) Values are based on magnitude differences directly measured by 
Daemgen et al. (2009). 



with a Kurucz atmosphere model to calculate a synthetic spec- 
trum. We then use STScI's pysynphot 2 to convolve the spectrum 
with transmission curves of all standard filter systems, taking 
into account the high-resolution Kepler response function 3 to 
obtain the throughput in Kepler 's broad spectral band. The re- 
sulting fractional contributions of third light, L3, are tabulated in 
Table 2. 



2.3. Fitting approach and error analysis 

To obtain adequate errors for the fit parameters and to determine 
the mutual dependence of the parameters, we explore the pa- 
rameter space by sampling from the posterior probability distri- 
bution using a Markov-Chain Monte-Carlo (MCMC) approach. 

Typically, the system parameters, especially the LDCs, are 
prone to substantial correlation effects. These correlations can 
easily render the sampling process inefficient. It is therefore ad- 
visable to use a suitably decorrelated set of LDCs as e.g. pro- 
posed by Pal (2008). However, due to the near-grazing system 
geometry in the case of TrES-2 we expect strong correlations 
between the whole set of parameters. We therefore chose to use 
a modification of the usual Metropolis-Hastings sampling al- 
gorithm, which is able to adapt to the strong correlation struc- 
ture. This modification, known as Adaptive Metropolis algo- 
rithm (AM; Haario et al. 2001), releases the strict Markov prop- 
erty of the sampling chains by updating the sampling parameters 
using a multivariate jump distribution whose covariance is tuned 
during sampling. This tuning is based on all previous samples 
of the chain, so that AM looses the Markov property. However, 
the algorithm can be shown to have the correct ergodic proper- 
ties, i.e. it approaches the correct posterior probability distribu- 
tion under very general assumptions (Haario et al. 2001; Vihola 
201 1). We checked that AM yields correct results for simulated 
data sets with parameters close to TrES-2, and found that this 
approach showed fast convergence and efficiency. 



3. Data preparation 

We retrieved the Kepler data of quarters Q0 to Q3 from the 
NASA Multimission Archive at STScI 5 . Specifically, we use 
data pertaining to data release 5 (Q0, Ql), data release 7 (Q2), 
and data release 4 (Q3) as provided by the Kepler Data Analysis 
Working Group 6 . Kepler observed TrES-2 in short cadence 
mode with a sampling of 58.85 s covering 229 days. For our 
analysis, we use the raw aperture photometry (SAP) flux and 
the corresponding error as provided in the FITS files. The raw 
data have been processed using Kepler s Photometric Analysis 
(PA) pipeline, which includes barycentric time correction, de- 
tection of Argabrightening (a transient increase of background 
flux possibly due to dust particles) and cosmic ray events, back- 
ground removal, and the computation of aperture corrected flux. 
Further, we removed invalid data points and all points marked by 
the SAP_QUALITY flag. The Kepler data provide time stamps 
in Barycentric Julian Dates (BJDutcX which we convert into 
Barycentric Dynamical Time (BJDtdb) accounting for 34 leap 
seconds elapsed since 1961. 

The data show discontinuities and exponential jumps due to 
instrumental duty cycles ("pointing tweaks", "safe mode recov- 
eries", and "earth points") and long-term trends ("focus drifts"), 
which must be corrected beforehand. First, we removed all tran- 
sits from the light curve including one half of the transit duration 
on either side and rejected 3<x outliers using a sliding median 
filter of window size 30 minutes. We then used the informa- 
tion provided in the data release notes to exclude data points 
affected by safe mode recoveries and earth points, which are 
difficult to correct otherwise. Subsequently, we divided the data 
into chunks covering undisturbed duty cycles. As there are two 
pointing tweaks occuring during quarter Q2, we subdivided the 
Q2 data at these tweaks into three chunks covering full duty cy- 
cles. 

To remove long-term trends due to focus drifts, we followed 
an approach similar to Kipping & Bakos (2011b) by using the 
discrete cosine transform (Makhoul 1980) to obtain a smoothed 
model of the continuum for each data chunk separately. The dis- 
crete cosine transform decomposes a signal x(n) into a linear 
combination of N cosine functions according to 



y(k) = J] x(n) 



cos 



xk(2n + 1) 
2N 



< k < N, 



(2) 



and provides an efficient way of extracting the low-frequency 
information within the signal. We applied the discrete cosine 
transform to each chunk of the light curve removing all but the 
first / low-frequency terms, where / is the rounded integer value 
of the length of the chunk divided by the orbital period. The 
thus cleaned transformed signal was inversely transformed to 
obtain a smoothed model of the continuum. Finally, we divided 
the complete light curve including the transits by the model 
and extracted the data in a box of half-width 0.06 days around 



2 http : // stsdas . stsci . edu/pysynphot/ 

3 http://keplergo.arc.nasa.gov/ 



CalibrationResponse . shtml 



4 http : //www .hs . uni-hamburg . de/DE/Ins/Per/Czesla/ 
PyA/PyA/ index . html 

5 http : //archive . stsci . edu/kepler/ 

6 http : //archive . stsci . edu/kepler/release_notes 
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Table 3. Publicly available data of ground-based photometric 
observations of TrES-2. The epoch refers to the first transit of 
TrES-2b observed by O'Donovan et al. (2006). 



Epoch 


Filter 


Reference 


13, 15, 34 


z' 


Holman et al. (2007) 


263 


near R 


Mislis & Schmitt (2009) 


312 


I 


Mislis & Schmitt (2009) 


395 


I 


Mislis etal. (2010) 


414 


l,y,b,v 


Mislis etal. (2010) 


421 


I 


Scuderi et al. (2010) 



each transit center. We further discard three transits, which occur 
close to the pointing tweaks, where the continuum flux shows 
strong distortions. 

To check that the results of our study do not depend on the 
details of the data reduction, we also extracted the un-detrended 
data directly around each transit center and divided the transit by 
a third-order polynomial fit to the local continuum (Welsh et al. 
2010). We find, that this procedure yields essentially the same 
results. 

The ground-based data covering a total of 8 planetary tran- 
sits between 2006 and 2009 are used as provided by the au- 
thors (see Table 3, for references). The time stamps are specified 
as Heliocentric Julian Dates (HJDtjtc) an d aie converted into 
BJDtdb using the webtool by Eastman et al. (2010). 

4. Fitting results 

Before fitting the Kepler data, we use all available data to obtain 
estimates of the period, P, and the time of transit minimum, t. As 
already pointed out by Southworth (2008), P and t are usually 
correlated with each other, but uncorrected with the remaining 
parameters. Therefore, the period can be determined quite accu- 
rately without interference of other parameters. Fitting all pub- 
licly available data we obtained 

T = 2453957.635486!°;™ BJDtdb and 

P = 2.4706 13402_ 0000000154 d, 

where to refers to the first transit of TrES-2b observed by 
O'Donovan et al. (2006) and errors correspond to 68.3 % HPD 
(highest probability density) intervals. The fractional error in P 
is approximately 10 5 times smaller than the fractional errors in 
a/R s and p, and can safely be ignored. We checked that includ- 
ing both parameters with Gaussian priors corresponding to their 
errors in the MCMC calculations yields negligible changes in the 
remaining parameters but slows down the MCMC algorithm. In 
the following, we therefore fix these parameters to their best-fit 
values. 

4.1. Four quarters of Kepler data 

In the following, we fit all 81 transits observed by Kepler 
during the first four quarters of observations. To explore the 
whole parameter space of solutions, we relax the constraints 
on the parameters as much as possible by imposing uniform 
prior probability distributions allowing large parameter varia- 
tions. By fixing the LDCs to their theoretically computed values 
we would impose strong a priori information on the fitting pro- 
cess. If the LDCs were chosen incorrectly, we would obtain well- 
constrained model parameters with comparably small credibility 
intervals, which, however, could be inconsistent. We therefore 



additionally fitted quadratic LDCs, assuming uniform priors on 
u\ and i<2, and performed an MCMC error analysis of the com- 
plete set of parameters. Our best-fit solution is determined by 
those parameter values that minimize the deviance after 10 6 iter- 
ations of the sampler. The errors are calculated from 68.3 % HPD 
intervals of the posterior probability distributions for the param- 
eters. Our results are summarized in the upper part of Table 4. 

Given their credibility intervals the parameters of our global 
fit are compatible with those by Southworth (201 1) and Kipping 
& Bakos (201 lb). We note a slight discrepancy with the results 
by Kipping & Bakos (2011b), who obtain a slightly larger ra- 
dius ratio, p, of 0. 1282. A small, possibly systematic discrepancy 
compared to the values of Kipping & Bakos (2011b) has also 
been noted by Southworth (201 1). Within their 2cr error intervals 
our LDCs are still consistent with the best-fit solution obtained 
by Kipping & Bakos (2011b), who found LDCs of u x = 0.52 
and m = 0.06. We note that the sum, u\ + m, of the LDCs is 
very well constrained in agreement with our considerations in 
Sec. 2.1. 

As discussed in Howarth (2011b) we cannot directly com- 
pare our best-fit LDCs to predictions based on the stellar atmo- 
sphere model. This is due to the fact that we are fitting pho- 
tometric data from a real star, using a simplified model of the 
stellar limb darkening. The high impact parameter of TrES-2b 
together with the fitting approach inevitably lead to consider- 
able deviations of the photometrically determined LDCs from 
model-atmosphere characterizations. Instead, Howarth (2011b) 
proposes to compare the results from the photometric analysis 
with the results obtained from synthetic photometry based on 
a given atmosphere model (termed "SPAM"). Figure 1 shows 
the result of this consideration. The lowest deviance combina- 
tion of LDCs after 10 6 MCMC samples is obviously inconsis- 
tent if compared directly to the predictions by Claret & Bloemen 
(2011) based on PHOENIX model atmospheres for TrES-2. 
However, the SPAM model based on the same model atmosphere 
is consistent with the best-fit LDCs. This strengthens our confi- 
dence in the reliability of the theoretically calculated LDCs. 

Note that the simplified prescription of limb darkening used 
in our model also introduces small, but non-negligible system- 
atic errors in the remaining light curve parameters due to the 
strong correlations between all parameters. In the following we 
therefore use the four-parametric, non-linear limb-darkening law 
proposed by Claret & Bloemen (201 1), holding the LDCs fixed 
at their theoretically expected values to obtain the remaining 
transit parameters. This, of course, substantially reduces the er- 
rors on our parameter estimates because the strong correlations 
mediated by the limb-darkening treatment diminuish. Our resuls 
after 10 6 iterations of the MCMC sampler are summarized in the 
lower part of Table 4. 

4.2. Secular change in inclination 

The values inferred from the fit to the Kepler data are the most 
robust estimates of the planetary parameters of TrES-2 currently 
available. In the following, we reanalyze all publicly available 
photometric data of TrES-2 examining each transit separately to 
search for secular variations of the orbital inclination. 

Secular inclination changes are empirically detectable as 
changes in transit duration, which is a function of a/R s , p and 
i (Seager & Mallen-Ornelas 2003). While each of these param- 
eters may itself exhibit secular changes resulting in changes in 
the transit duration (Kipping 2010), we focus in the following on 
changes in the orbital inclination, i, and the time of transit min- 
imum, t, of each transit. We thereby assume that neither the ra- 
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Table 4. Lowest-deviance parameters for the Kepler data after 
10 6 MCMC samples compared to the results with fixed limb- 
darkening coefficients. Both fits assume a third light contribution 
according to Table 2. 



Parameter 


Value 


Error (68.3 % HPD) 


in 


83.874 


[ 83.806,83.910] 


p 


0.1257 


[0.1240,0.1265 ] 


a/R s 


7.903 


[ 7.835 , 7.937 ] 


"i 


0.19 


[-0.10,0.36] 


"2 


0.43 


[ 0.25 , 0.75 ] 




0.626 


[0.605, 0.648] 


hi 


83.9788 


[ 83.9740 , 83.9837 ] 


p 


0.12759 


[0.12755 ,0.12762] 


a/R s 


7.9899 


[ 7.9856 , 7.9943 ] 




fixed, see Table 2 




-0.2 - 



-0.4 1 1 1 1 L - 1 1 1 1 

-0.4 -0.2 0.2 0.4 0.6 0.8 1 

Coefficient Uj 

Fig. 1. Correlation between linear and quadratic limb-darkening 
coefficient resulting from 10 6 MCMC runs. The lowest deviance 
combination of parameters (circle) is significantly offset from 
the Claret & Bloemen (201 1) model prediction (dashed red) but 
consistent with the result obtained from fitting synthetic photom- 
etry based on the same atmosphere model (SPAM, green). 

Table 5. Prior information used for the MCMC calculations in 
Sec. 4.2. 



Quantity 


Prior 


Source 


a/R, 


7.9899 ± 0.0043 


Table 4 


P 


0.127585 ± 0.000035 


Table 4 


T 


uniform 




i 


uniform 




fill, Cl2, Cl3, an 


fixed 


Table 2 




fixed 


Table 2 


e 


(fixed) 





dius ratio nor the planetary orbit change in size over the course of 
the observations and that the orbit is circular. However, we still 
must account for our incomplete knowledge of the true value of 
a/R s and p. We therefore leave both parameters free during the 
MCMC calculations imposing Gaussian priors according to the 
errors obtained from our global fit. Depending on the instrumen- 
tal setup of the observations (cf. Table 3) we take into account 
the expected third light contribution and the LDCs according to 
Table 2. For each transit we obtain the posterior distributions 
for both t and ; using 10 4 iterations of the MCMC sampler. We 
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Fig. 2. Inclination versus transit epoch for TrES-2 based on pub- 
licly available ground-based data and the Kepler observations. 
Errors are derived from 95.5 % HPD intervals corresponding to 
2cr errors based on 10 4 MCMC iterations for each transit. 



summarize the prior information in Table 5. The errors again 
correspond to 68.3 % HPD intervals. 

Concerning the Kepler data, Kipping & Bakos (2011b) did 
not find evidence for parameter changes within the first two quar- 
ters. However, the public release of the last two quarters, i.e. 
quarters Q2 and Q3, provided a substantial enlargement of the 
data set. We therefore reanalyze the complete Q0 to Q3 set of 
transits in detail searching for variations of the inclination or 
transit timing variations (TTVs). 

Our results are shown in Fig. 2. The superior quality of the 
Kepler data manifests itself in very small credibility intervals 
for these observations. While the data by Holman et al. (2007) 
are fully consistent with the inclination obtained by Kepler, we 
find systematically lower inclinations for the data sets of Mislis 
& Schmitt (2009) and Mislis et al. (2010). In fact, ignoring the 
Kepler data, the secular inclination change detected by Mislis & 
Schmitt (2009) is clearly visible. 

Note, however, that our common fitting approach based on 
the results from Kepler and including the third light contribu- 
tion yields best-fit values for the inclinations derived from the 
ground-based measurements that are closer to the Kepler values 
than those originally obtained by Mislis & Schmitt (2009) and 
Mislis et al. (2010). Furthermore, our new errors are quite large 
so that the Kepler value lies only just outside of our 3cr error. 

We also checked, whether it is possible to completely recon- 
cile the Kepler observations with the results by Mislis & Schmitt 
(2009) and Mislis et al. (2010). Accounting for the possibility 
of normalization problems during data reduction we introduce 
an additional normalization constant to our model. Including 
this constant as another parameter of our MCMC marginally en- 
larges the uncertainties on the inclinations. We thus do not ex- 
pect the deviation to stem from simple normalization problems. 

Another cause for erroneous results would result from in- 
correctly taking into account the contribution of the contaminat- 
ing third light. To fully reconcile the inclinations by Mislis & 
Schmitt (2009) and Mislis et al. (2010) with the best-fit inclina- 
tion from the Kepler data would require contributions of a third 
light of 10 to 20 % of the total flux, which can be excluded from 
the observations by Daemgen et al. (2009). However, assuming a 
uniform prior for the third light contribution and demanding it to 
be smaller than 10%, the credibility intervals become substan- 
tially larger including the Kepler result. We thus conclude, that 
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problems with the flux calibration possibly caused by the refer- 
ence star are the most likely cause for the observed deviations. 

In either case, the Kepler data rule out a decrease in inclina- 
tion of the size expected by Mislis et al. (2010). The exact reason 
for their systematic deviation remains unclear. The Holman et al. 
(2007) data are consistent with the first four quarters of Kepler 
data. Note, however, that despite the long available time baseline 
the errors of the Holman et al. (2007) data points are too large 
to draw significant conclusions on the presence of other, smaller 
linear trends. 

4.3. Variability in the Kepler data 

We now focus on the Kepler data searching for small trends 
or variations, which could still be consistent with the previous 
ground-based observations. Fig. 3 shows an enlarged view of 
the Kepler data; the errors are derived from the 68.3 % HPD in- 
tervals, so they are equivalent to 1 cr errors. 

We found a slight linear trend in the Kepler data. 
An MCMC analysis yields a slightly non-zero slope of 
(0.8 + 0.2) x 10~ 4 "/cycle. This would correspond to a secular 
inclination change nearly two orders of magnitude smaller than 
the value by Mislis et al. (2010). Extrapolating this trend back 
to the epoch of transit observations obtained by Holman et al. 
(2007), we would expect an inclination of 83.94° ±0.01° consis- 
tent with the measurement errors. 

In order to quantify the improvement of the different fits, we 
calculated the statistic 

f - — (3) 

xiln 

and compared its value to an F-distribution with v a - Vb and v\, 
degrees of freedom (e.g., Rawlings et al. 1998). We obtain a p- 
value of 0.004 indicating that the linear trend provides a better 
description compared to the constant at the 2.8cr level, which 
we consider as indicative, but not significant evidence for the 
presence of such a trend. 

After subtracting the linear trend from the data, we computed 
an error-weighted Lomb-Scargle periodogram (Zechmeister & 
Kiirster 2009), which displays a single signal at a period of 45 
cycles or 111 days with a false-alarm probability below 10%. 
Fig. 3 shows our best-fit sinusoidal model including a linear 
trend. An F-test comparing the model to the constant model 
yields a p-value of 0.015 (or 2.4cr); the test thus indicates the 
presence of an additional sinusoidal modulation superimposed 
on the linear trend, which however remains insignificant. 

We can also compare the three models using the Bayesian 
information criterion, BIC - x 1 +k\s\N (Schwarz 1978), which 
penalizes the number k of model parameters given N - 2,1 data 
points. We obtain BICs of 136.3, 127.6, and 117.0 for the con- 
stant, the linear trend, and the sinusoid with linear trend, respec- 
tively. 

Despite being close to formal significance the signal is eas- 
ily obliterated by the noise and could be introduced by some 
subtleties in our analysis. Assuming the sinusoidal model with 
linear trend to be true, we averaged two sets of 10 transits per- 
taining to low (around epoch 409, compare Fig. 3) and high 
(around epoch 470) inclinations, respectively. Subtracting both 
templates, their difference should exhibit a characteristic signa- 
ture, which we model by subtracting two model calculations for 
i'z, = 83.97° and in = 83.99°. We show the averaged templates 
together with their difference and the model in Fig. 4. While 
during ingress the flux difference shows a decrease following 



the difference model, the signature during egress is obliterated 
by the noise. We checked that averaging transits pertaining to 
epochs close to the zero-crossing points of the sinusoidal model 
yields no substantial deviation from a constant. 

We repeated our analysis using simulated control data to 
check, whether the signal is introduced by our data reduction 
or by "phasing effects" due to the Kepler cadence (Kipping & 
Bakos 201 la). However, we could not detect similar signals ex- 
ceeding the noise. We further changed our data normalization 
strategy (see Sec. 3), but found the variation of inclination to be 
stable against these changes. Analyzing similar Kepler targets 
we found spurious periodicities in some cases, which, however, 
remained insignificant in the periodogram (at false-alarm prob- 
abilities above 10 %). In summary, there seem to be some indi- 
cations for a short-term inclination variation and an underlying 
secular change in the Kepler data of TrES-2. However, the de- 
tected signal is merely an indication of inclination changes in the 
Kepler data, because in the baseline of space-based observations 
is still too short to draw significant conclusions on the stability 
or variability in the inclination (compare Fig. 4). 

Turning to the times of transit minimum of the Kepler tran- 
sits, Fig. 5 shows that the transit timing is consistent with being 
constant over the course of the Kepler data. An F-test compar- 
ing the constant and a linear model yields a p-value of 0.92; thus, 
there is no evidence for a linear trend. In the periodogram there 
is no signal at false-alarm probabilities below 10 %. We thus do 
not see any evidence for TTVs in the Kepler data. 

The apparent secular change in inclination refers to the ori- 
entation of the orbital plane of TrES-2b relative to the observer's 
tangential plane and, assuming a circular orbit, can be inter- 
preted as nodal regression at fixed orbit inclination. While grav- 
itational perturbation due to an oblate host star is unlikely given 
that TrES-2A is a slow rotator very similar to the Sun (Mislis 
et al. 2010), additional bodies in the stellar system could offer a 
viable explanation for secular changes in the orbital parameters. 
If we assume our linear trend in inclination to be physical, we 
can estimate the period of a perturber of specified mass in a cir- 
cular, coplanar orbit using Eq. 13 in Mislis et al. (2010). We find 
that a perturbing planet of Jovian mass would be expected at or- 
bital periods between 100 and 300 days. However, a perturbing 
second planet capable of causing nodal precession of the transit- 
ing planet would also be expected to induce short-term periodic 
variations of its orbital elements (Holman & Murray 2005; Agol 
et al. 2005). To calculate the amplitude of such TTVs, we use 
an approach based on analytic perturbation theory developed by 
Nesvorny & Morbidelli (2008). Outside possible mean-motion 
resonances additional Jovian mass bodies with periods larger 
than 100 days would cause a TTV amplitude below 1 s if on a 
circular orbit. Given the standard deviation of the times of transit 
minimum of 4.8 s, such a TTV amplitude cannot be ruled out. 
Thus, the existence of an external perturber in the system still 
provides an entirely plausible scenario. 

5. Summary and Conclusion 

We presented an analysis of all publicly available photometric 
data obtained during planetary transits of TrES-2b including four 
quarters of data obtained with the Kepler space telescope. The 
transits were fitted in a common scheme using theoretically cal- 
culated LDCs and taking into account the contribution of a third 
light. 

Fitting quarters Q0 to Q3 of the Kepler data, we were able 
to refine the system parameters determined by Kipping & Bakos 
(201 lb) and Southworth (201 1). Our results are compatible with 
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Fig. 3. Inclination versus transit epoch for the Kepler observa- 
tions. Errors are computed from the 68.3 % HPD intervals based 
on 10 4 MCMC iterations for each transit. The best-fitting lin- 
ear (dotted/green) and sinusoidal models including a linear trend 
(dashed/red) are overplotted. 
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Fig. 4. Upper panel: Averaged templates of two sets of 10 tran- 
sits with high (H, solid red) and low (L, dashed green) incli- 
nations, respectively, as expected from the tentative sinusoidal 
variation of the inclinations obtained by Kepler. Lower panel: 
The difference between the two averaged transits. The expected 
difference for two models pertaining to ; = 83.97° and 83.99° is 
shown for reference (solid red). 



these previous studies, our errors being slightly smaller. Its high 
impact parameter makes the TrES-2-system sensitive to details 
of the adopted stellar limb-darkening model because only the 
covered regions close to the stellar limb are used to extrapolate 
the limb darkening for the whole stellar disk. Small differences 
between limb-darkening prescriptions therefore yield significant 
differences in the resulting best fits. We find that the theoretically 
calculated LDCs by Claret & Bloemen (201 1) provide a correct 
description of the observed limb darkening if compared with the 
data via the SPAM approach as advocated by Howarth (201 lb). 

The data by Mislis & Schmitt (2009) and Mislis et al. (2010) 
are found to be inconsistent with the remaing data, especially 
the Kepler data. The discrepany can neither be attributed to er- 
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Fig. 5. Transit timing variation (TTV) versus transit epoch for 
TrES-2b using the best-fit ephemerides. Errors are computed 
from the 68.3 % HPD intervals based on 10 4 MCMC iterations 
for each transit. 



rors in the theoretical LDCs nor to the contribution of third light, 
which the original studies did not account for. The most proba- 
ble explanations would be subtleties in the data reduction, per- 
haps due to problems with the reference star, or the presence of 
time-correlated noise yielding underestimated parameter errors 
(Carter & Winn 2009). 

Excluding the Mislis & Schmitt (2009) and Mislis et al. 
(2010) data sets, all observations can be reconciled in a consis- 
tent picture. The extension of the available Kepler data by two 
additional quarters reveal some hints for systematic variations in 
the inclination, the marginally significant slope would be consis- 
tent with the data by Holman et al. (2007). However, given the 
available baseline of space-based observations we consider the 
inclination change as indicative but not significant. Further data, 
which already have been obtained by the Kepler space telescope, 
will settle this issue in the near future. 

Our study shows that space-based and ground-based pho- 
tometry of TrES-2 can be analyzed in a common scheme yield- 
ing consistent results. While space-based photometry offers the 
possibility to constrain all system parameters at once, with the 
eventual decommissioning of Kepler it will again be up to 
ground-based measurements to search for long-term changes in 
the orbital parameters of TrES-2. 
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Appendix A: Results from individual fits 

In Table A.l, we show the transit parameters of TrES-2b from 
individual fits of ground-based and Kepler data. 
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Table A.l. Transit parameters for TrES-2b from individual fits of ground-based and Kepler data. Epoch refers to the first transit of TrES-2b 
observed by O'Donovan et al. (2006). For the MCMC analysis, we fix the limb-darkening coefficients corresponding to the filter band (cf. Table 2) 
and use prior information as described in Table 5. 
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1 97,8O+0000° 37 
1Z./JOU 000031 


419 


499? 8??5?7 +0 000042 
t77z.Oiz.ji / _ .000052 


83 


O040+O.OO54 
-0.0062 


7 


080^+0.0048 
,7 ° -0.003! 


0. 


1 ?7S7fS+0000 036 
—0.000035 


420 


4995 293008 + °° 00044 

t77J.z.7JUUO n .000046 


83 


Q7 if. +0.0051 
■ y 1 1 -0.0057 


7 


989O+0- 0041 
.yoyv_ QQQ40 


0. 


1 97S78 +0-000 035 
—0.000032 


421 


4997 763604+ 0000045 

^yyi.i «-'«^' t _o.000045 


83 


Oonc+0.0051 
■ yoKJJ -0.0062 


7 


O8QS+0.0039 
■' o,J -0.0043 


0. 


1 97S78 +0-00 0038 
—0.000032 


424 


5005 175613+ 0000049 

juuj.i /jui J 0.000044 


83 


0741+0.0063 
■ y '^ L -0.0049 


7 


99OO+OOO 41 
.yy\j\j _ 00041 


0. 


1 97S 89 +0 00 0033 
iz,/ J 0:7 -0.000036 


425 


5007 6461 98 +0 000045 

juu/.utui?o 0.000049 


83 


Q81 O+0.0053 
■ yM8 -0.0056 


7 


9884+00043 

■ ;?OC ' H -0.0041 


0. 


1 97,77+0.000037 
/ -0.000032 


426 


5010.116754!°;™ 


83 


071 c +0.0055 
■ ylLJ -0.0056 


7 


99 04+0.0043 
yyu -0.0039 


0. 


1 97<C8^+0-0000 40 
iz./joj_0.000029 


427 


5012.587380!";™ 


83 


9809+0.0056 

■yo\jy_ QQQ59 


7 


Q8QQ+0.0045 
■ :?o: ' : '-0.0041 


0. 


1 97, 8, +0-000030 
—0,000038 


430 


5019 999187+ 0000042 

-0.000052 


83 


9850 +0.0055 
■ 7OJU -0.0060 


7 


989S +0.0048 
.?o? J_ 0004 5 


0. 


1 77S 89 +0 00 0032 
iz./jo5'_o.ooo036 


431 


5022 469822+ 000051 

juzz.+u?ozz_ .000044 


83 


97 Oo+o.oo52 

,7/u -0.0058 


7 


9901 +0.0039 
,: ' : ' L ' 1 -0.0044 


0. 


1 97S 89 +0 00 0037 
iz./jo5'_0.000033 


432 


S0?4 940477+0 0000 50 
juzt.7W+ / / _ .000044 


83 


077 1 +0.0051 
■ y 1 ' 1 -0.0067 


7 


Q8Q9 +0.0043 
.^O^i_ooo46 


0. 


1 77S 81 +0-000040 
000030 


433 


5027 411097+ 0000046 

juz/.tiiU7( -0.000047 


83 


O 862 +0.0050 


7 


9897+00041 
■ yoyl -0.0041 


0. 


1 97 ,79+0 000° 35 
1Z./J '^-0.000036 


434 


5029 881709+00000 43 

JUZ, - M1 /u ^-0.000048 


83 


9768 +0.0054 
■ ?/uo -0.0055 


7 


0007+O.OO4O 
-0.0043 


0. 


1 97,77+0.000034 
/ -0.000037 


435 


SfVo9 a<99 1 9 +0.000047 
^-0.000045 


83 


0744+00057 
y -0.0058 


7 


0809+0.0041 
,:7o: ' : '-0.0044 


0. 


1 97,8O+0000° 38 
li,JO -0 .000032 


436 


5034.822785!°;™ 


83 


98,n+0°0 52 
.?oju_ 00058 


7 


99Q6 +0.0039 
'^""-0.0045 


0. 


1 97<C87 +0-000° 32 
000035 


437 


5037.293526!";™ 


83 


9665 +0.0052 
■ 7UUJ -O0060 


7 


Q8QQ+0.0041 
.yoyj_ Q(m2 


0. 


1 97S77+0 000 034 
' -0.000036 


438 


5039 764149+ 000046 

JUJ?.)OTit7 -0.000049 


83 


0749+0.0057 
■ y -0.0056 


7 


989O+0 0044 
■ ; ' o; ' u -0.0041 


0. 


1 97S 87 +0 00 0031 
1Z/JOZ -0 .000038 


439 


5042 234769+ 000049 

jutz..z. / u^o.000047 


83 


08 1 q+0.0054 
,?01 -0.0063 


7 


QQ 16+0.0042 
,?: ' 1U -0.0046 


0. 


127606+ 0000035 

iz.;ouo_o .000036 


440 


5044 705462+ 000049 

jutt. / ujtu^ o.000046 


83 


Q701 +0.0061 
-0.0054 


7 


9897+00041 

-0.0044 


0. 


1 97S8O+0000 037 
ow_o 000034 


441 


5047 175958+ 0000043 

jut / . 1 / j^jo_ .000047 


83 


9781 +00059 
-0.0054 


7 


9889+ - 0038 
.yooy_ QQQ41 


0. 


1 ?7,8V0-000 034 
lz.1 JO J -0.000035 
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Epoch r i 


a/R s 


P 


Reference 


[BJDtdb - 2, 450, 000 d] [deg] 









442 


5049 


f.Af.z 1 o +0.000046 
.utuj±o_ 0000050 


83 


qf>qi +0-0057 

■ y " yL -0.0050 


7 


0007+00039 
-0.0042 


0. 


1 77CQ1 +0.000038 
LZ -'~> yL -0.000032 


443 


5052 


1 1 71Q4+0.000045 


83 


Q774+0.0063 
■ yl -0.0047 


7 


0007+00035 


0. 


1 77<88+°-000031 
iz.; ~>°°_o.000038 


444 


5054 


S 8774.7 +0-000051 
/H/ -0.000046 


83 


077C+O.OO58 
• yl /J -0.0057 


7 


Q88S+^-^ 44 
-0.0042 


0. 


1 ?7S 81 +0-000040 
li/JO -0.000030 


445 


5057 


r\ c o o £L o +0. 000046 
.ujojuo_ 0000049 


83 


O/;7^+0.0055 
■ yu,J -0.0054 


7 


Q8QO+ - 0044 
■ ?o: ' u -0.0039 


0. 


1 ?7S81 +0.000037 
-0.000034 


446 


5059 


S?Q1 <<+0.000044 
' J ^ 1J -0.000049 


83 


Q707+O.OO56 
■ y,c " -0.0060 


7 


0800+00042 
yoyyJ -0.0045 


0. 


1 77S80+^-^^ 3 ^ 
IZ ' J -0 .000033 


447 


5061 


QQQAC 1 +0.000048 
■ yyyK >oi _ 000045 


83 


Q76Q+00054 

-0.0057 


7 


Q8QS+ 0043 
-0.0042 


0. 


1 77<;7/;+0-000036 
iz,/J/u -O.O0O035 


449 


5066 


Q/l OB Q1 +0.000048 

,940882_ 0000045 


83 


Q786+00050 
■^ /ou -0.0060 


7 


QQ07 +0.0043 
■" u -0.0041 


0. 


1 ?7S8q +( ' () "^ 34 

uiJO -0.IK10036 


450 


5069 


41 i4Sfi +0000046 


83 


Q70C +0.0050 
■ :?/ - : ' -0.0058 


7 


9896+0.0042 

■^°^ u -0.0041 


0. 


1 77<8^+ ( '- ( '"^ 33 
iz.; joj_ .000036 


451 


5071 


8 8700? +0-000045 


83 


qs 14+0.0051 
-0.0061 


7 


QQOO+ 0048 

■"""-0.0038 


0. 


1 97SR7+0-000033 
iz. / jo / —0.000035 


452 


5074 


oc97/;£:+0. 000044 
JJZ/UU -0.()00048 


83 


9856 +00058 

■ 7OJU -O0054 


7 


qq 1 o +0 0040 

■ yyl -0.0046 


0. 


1 ?7Sqi +0.000034 
U/J:,1 -0.000033 


453 


5076 


o 9 o o A 7 + ■ 000047 
■ oz - JJ ^ '-0.000048 


83 


9805 + 00058 

,7OUJ -O0050 


7 


qoq 6 +0.0045 
yoy -0.0035 


0. 


1 97V7B+0.000037 
iz./j/o_o.ooo032 


454 


5079 


9QOQ7A+0. 000049 
.z,^j^/w_ 0000043 


83 


n 7 g 6 +0.0051 
y '° -0.0058 


7 


QQO4+00042 
.?7Ut_ .0042 


0. 


1 ?7S8 1 +0.000039 
000029 


455 


5081 


■ A AAAA A 7 

764470+ 0000043 
,/utt/u Q 000054 


83 


O7SQ+00051 
•^'^ -0.0059 


7 


1 /a f\n 1 1 

080? 

■'°' z -0.0O38 


0. 


1 A AAAAT C 

1 ?7S8 1 +0.000035 
iz. / jo 1 0.000035 


456 


5084 


7-3 C74C+0. 000045 
• iJJ ^ J -().0()0048 


83 


qqfi5+0.0058 
,: " OJ -O0048 


7 


qq07 +00037 
■ yyK)l -0.0043 


0. 


1 77 <Q4+0. 000036 
li,J: ' -0.000032 


457 


5086 


70^707 +0.000048 
''"-"^ -0.000044 


83 


Q8 7 7 +0.0053 
.yojz._ 0Mm 


7 


QQO? +0.0045 
-0.0040 


0. 


1 ?7^8 1 +0.000035 
li/JO -0.000034 


460 


5094 


1 1 754A+0.000046 
.11 '-'' + "_o.000047 


83 


q744+0.0065 
■^'^ -0.0046 


7 


00 1 4+0.0041 
yyL -0.0043 


0. 


1 77 <Q4+0. 000037 
li,J "-0.000O33 


461 


5096 


c +0.000042 
.joo/,ro_ 0(|(|()(|49 


83 


Q7Q4+0.0057 
■ 7/7 -0.0052 


7 


QQO4+00042 
•" u -0.0039 


0. 


1 77S8S+^-^^ 32 
1 z.1 JO j_ Q ooo038 


462 


5099 


C ,8744 +0000049 
' UJO '^ -0.000046 


83 


QB7<; +0.0059 
-0.0057 


7 


0008 +00042 
■"""-0.0048 


0. 


1 77SQ7+0-000037 
iz./j^j_ ()(K)0()34 


463 


5101 


c 9 q o n n +0. 00005 1 
.jz?j/u () 000044 


83 


984S+ 00058 

.=-otj_o.0051 


7 


QQ07 +0.0043 
.7?UJ_ .0043 


0. 


1 77SQQ+0000037 
iz./ J77_| ) . 0000 32 


464 


5103 


qqqqqO+O- 000048 
,7777 -O.IM)M5 


83 


Q7Q7 +0.0050 
y,y -0.0058 


7 


qsqs+ooo 43 

■ 7OJJ -0.0041 


0. 


1 ?7 Sq? +0.000037 
li/J7i -O.OO0O33 


465 


5106 


47O7OfS +0000049 
-t/u/uu Q .000043 


83 


qo<<+0.0049 
,,OJJ -0.0061 


7 


q8qq+ 0044 

■ :?o;?: '-0.0041 


0. 


1 ?7S8fi +( ' ( '" ( ' ( ' 35 
li ' JOU -0.000035 


466 


5108 


q41 77S+0. 000046 
-0.000044 


83 


004^+0.0056 
.^?tj_ 061 


7 


q8q7+0°048 
■^ 0: " -0.0037 


0. 


1 77C8 1 +0.000038 
000030 


467 


5111 


41 1 BOB +0.000042 


83 


q7qfi+0.005l 
■ 7/7o -0.0059 


7 


9896+0.0041 

■^°^ u -0.0044 


0. 


1 77<8^+°-000035 
iz./ joj_ .000037 


469 


5116 


o c-i 078 +0.000042 
.jjju/o_ 0(K)004g 


83 


no 1 a +0.0062 
-0.0049 


7 


QQO7+0 0046 
yyy " -0.0039 


0. 


1 ?7Sqi +0000035 
"'■"-0.000033 


470 


5118 


0700^7+0.000046 
.oz,jouz,„ ()00049 


84 


0001+ 00058 
-0.0051 


7 


QQ77+0.0047 
.??z.z._ 00038 


0. 


1 77«;qq+0000033 
iz/ J^^_ .000038 


471 


5121 


7Q47QI +0.000047 
■"^-0.(100046 


83 


Q077+O.OO53 
yc " -0.0056 


7 


QSOQ+O- 0040 
^-0.0046 


0. 


1 77<;qo+ ( ' ( '0^ 32 
li,J -0 .000036 


473 


5126 


ttccoi +0.000044 
' ZJJJOJ -0.000047 


83 


q8 8 O+0.0052 
■ y888 -0.0061 


7 


Q8Q7+0.0043 
.?o?j_ 0040 


0. 


1 ?7S8fi +( ' ( '"^ 35 
iz/ jou _ 000037 


474 


5128 


■ a aaaa a a 

7061 q? +0.000046 

./uui?z._ 0000046 


83 


qq41 +0 0051 

yy ^ -0.0058 


7 


1 A A A A A 

qq 1 o+o- 0044 

,771U -0.0036 


0. 


1 77<8^i+0-000037 
liM -0 .000033 


475 


5131 


1 768^4+0 000042 
.1 /UOJt_ Q 000047 


83 


Q877+O.OO6I 
' 7 -0.0045 


7 


QQ 1 1 +0.0040 
yy 11 -0.0039 


0. 


1 77<Q7+0.000034 
"'-^ 7 -0.000035 


476 


5133 


6474? ^ +0 000042 


83 


QQ7 1 +0.0052 
yyJL -0.0058 


7 


QQOQ+OOO 41 
■ :7:7U: '-0.0041 


0. 


1 77<Q7+0000033 
1Z/J,J -0.000035 


477 


5136 


1 1 BHBA+0. 000043 

•118086_ ()()00052 


83 


9800 +00055 

,?OUU -0.0059 


7 


0007+00042 
■ 77U/ -0.0046 


0. 


1 ?7SQ? +0-000029 
lz ' J,i -0.000043 


478 


5138 


^ QQ&HH +0.000048 
.5886//_ .000047 


83 


9806 + 00055 

.^ouu_o o60 


7 


OQ78 +0.0040 
■"^°-0.0049 


0. 


1 ?7607. +0 000035 
iz - /UUJ -0.000035 


479 


5141 


nsq?s4+°o° 0044 

•"•"^ -0.OOOM9 


83 


9840+ 00060 

.^ot -0.0050 


7 


QQOS+ 0041 
■ 7WJ -0.0043 


0. 


1 77«;88+ ( ' ( '0^ 38 
iz/joo _o.ooo031 


480 


5143 


c 7 q 7 + . 000045 
-jzj'oo/ _ .000043 


83 


Q7Q7+O.OO6I 
■ y,y ' -0.0049 


7 


QQCQ+0.0043 
■ y88 8 -0.0041 


0. 


1 77^77 +0.000036 
iz/j/ ' -0.000034 


481 


5146 


■ A AAAA \ ~> 

000644 +0(K)0()42 

• wuu -0.000046 


83 


Q<;<q+0.0053 
■ 7u -> -0.0059 


7 


1 A A/110 

080^+00038 


0. 


1 A AAA AO 7 

1 ?7S8O + 0'0O0033 
liMO -0.O0OO36 


482 


5148 


47 1 1 7 7. +0.000046 
.t/iijj_ 0()00046 


83 


08 1 7+0.0052 
■^ oi -0.0057 


7 


q8q8 +0 0043 


0. 


1 ? 7 S R4+0.000036 
li,J -OOOM33 


483 


5150 


q41fiqi +0 000043 
■ y ^ L " yl _ ,000047 


83 


q7qq+00()54 
•^'^ -0.0055 


7 


q8qq +00()42 

,:,o: ' : '-0.0039 


0. 


1 ?7SQ1 +0-000034 
U ' J '-0.000037 


484 


5153 


41 77 B^: +0.000044 
.tiz.z.oj_ 000045 


83 


Q7CC +0.0049 
■ y 1 J -0.0062 


7 


QQO? +0.0040 
■" u -0.0042 


0. 


1 ?7S8^ + ^-^^ 3 ^ 
iZ,/JOJ -0.000036 


486 


5158 


<on+0. 000045 
000049 


83 




7 


0000 +00043 
^'"'-O.MMO 


0. 


1 77SQ7+0-000034 
iz/J5»z_ 0000034 


487 


5160 


8?40q6+o° 00046 
.oz.tu^u_ .000043 


83 


07 +0.0054 
ylJJ -0.0062 


7 


0000 +00047 
.7?u?_ 004 o 


0. 


1 77SQ7+0000035 
li,J7i -0O00O35 


488 


5163 


. . _ 1 A AAAA 1 1 

?q47 1 7+0.000041 
' i,H ' JJ -0.1M0051 


83 


q746+°oo 56 

■ 7 '^"-0.0054 


7 


„ „ , A AM ^ O 

0007 +00042 
■ ;?: ' U/ -0.0042 


0. 


, l\ AAAA"} T 

1 77<87 

000037 


489 


5165 


765446+ 0000051 

. / u-'' + ' +< J_o.000044 


83 


Q87A+0.0058 
,70J -0.0053 


7 


QOQ7 +0.0047 
yoy 1 -0.0036 


0. 


1 77S 87 + - 000032 
iz,/JO/ -0.000035 


490 


5168 


236001 +0000042 

—0.000047 


83 


q 1 c ( : +0.0056 
,7/ JU -0.0058 


7 


qqi 7+0.0041 
■ ;?: ' i - : '-0.0045 


0. 


1 77<CQ7+0.000031 
iz/ j?j_ 0000036 


491 


5170 


706625 +0 - 000045 

•'" u " iJ -0.000043 


83 


Q7QO+0.0048 
■ 7 ' 7U -0.0061 


7 


QQOS+0 0046 
,77 " J -0.0038 


0. 


1 ?7SqO +0000036 
1Z ' J?U -0.000032 


492 


5173 


1 771 BQ+0.000049 
.1 / / i °5'_o.000043 


83 


go 1 4 +0.0057 
,7ol ^-0.0052 


7 


QQO7+0.0039 
■ :7:7U/ -0.0044 


0. 


1 ?7<;q8 +0000034 
iz/j^o_ 000003 5 


493 


5175 


^7077+0.000043 
-0.000051 


83 


967 o+0.0050 
yul °-0.0059 


7 


QQ0 1+0.0042 
■ yyKJL -Q.0Vh9 


0. 


1 ?7S81 +0.000033 
iz/ jo 1 0.000036 


494 


5178 


1 1 64^0 +0.000039 


83 


078 1 +0-0060 
■ y,OL -0.0055 


7 


QQO7+0.0048 
■ yy w -0.0043 


0. 


1 77sq8+00° 0034 
iz/j^o_ .ooo036 


495 


5180 


58q?04 +0000040 
J ° 7 -0.000053 


83 


Q74.7+0.0048 
■ ?/ ^' -0.0057 


7 


Q8Q7 +0.0043 
.?o?z._ 00037 


0. 


1 ?7S8? +0 - 000035 
U ' JOZ -0.0(I0036 



